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EQUATORIAL  SPREAD  F: 

“IN  SITU”  MEASUREMENTS  OF  ELECTRON  DENSITY 
TEMPERATURE  AND  DENSITY  FLUCTUATION  POWER  SPECTRA 

I.  INTRODUCTION 

The  DNA/PLUMEX  rocket  payloads  launched  into  the  equatorial 
ionosphere  during  the  July  1979  campaign  carried  a  plasma 
diagnostics  complement  that  included  a  quadrupole  ion  mass 
spectrometer,  a  pair  of  pulsed  plasma  probes,  vector  electric 
field  sensors  and  a  four-frequency  beacon.  The  pair  of 
pulsed  plasma  probes  not  only  provided  simultaneous  measure¬ 
ments  of  electron  density  Ng,  temperature  Tg ,  and  density 
fluctuation  power  spectra  P  ( k ) ,  but  also  provided  the 
capability  for  a  running  measurement  of  relative  variations 

in  mean  ion  mass  <M.>. 

1 

The  measurements  of  N  and  T  form  the  basic  information 

e  e 

on  the  laminar  condition  of  the  ionosphere,  allowing  for  the 
determination  of  the  ionospheric  plasma  response  to  varying 
geophysical  conditions  (solar  and  magnetic  activity,  winds, 
gravity  waves,  etc.)  and  the  detection  of  triggering  mecha¬ 
nisms  (e.g.,  steep  density  gradients  in  Ng)  for  ionospheric 
irregular  it  ies . 

The  measurements  of  <5N  (— •  P  (k))  yield  important  test 

e  n 

information  for  signal  channel  models  as  well  as  candidate 
instability  mechanisms  (e.g.,  collisional  drift  modes,  E  X  B 
and  Rayleigh-Taylor)  which  might  be  active  in  the  ionospheric 
p 1 a  sma . 

In  this  paper  we  describe  the  experimental  technique, 

payload  configuration,  launch  scenerios  and  present  additional 

1  2 

experimental  results  which  complement  companion  papers  ’  . 

Manuscript  submitted  June  13,  1980 
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II.  TECHNICAL  OVERVIEW 

3 

Introduction.  The  pulsed  plasma  probe  technique  (P  is 

the  designated  acronym)  is  a  unique  diagnostic  tool  capable 

of  high  spatial  and  temporal  resolution  of  plasma  parameters. 

The  instrument  is  a  Langmu ir- t y pe  probe  using  a  special 

electronic  procedure  for  generating  the  current-voltage 
3  4 

characteristic  ’  .  The  result  is  greatly  improved  reliability 

and  expanded  versatility  in  Langmuir  probe  measurements.  As 

3 

a  diagnostic  tool,  the  P  technique  reduces  commonly  found 

distortions  in  derived  electron  densities  and  energy  distribution 

functions.  A  unique  feature  of  the  technique  is  its  ability 

to  measure  simultaneously  the  electron  temperature,  density, 

and  the  density  fluctuation  power  spectrum.  Successful 

3 

applications  of  the  P  technique  include  not  only  rocket  but 
also  satellite  5,6  and  laboratory  beam-plasma  studies7  of 
turbulent  charged-particle  environments. 

Figure  1  shows  two  types  of  probe  operation.  Figure 
1(A)  depicts  a  linear  sawtooth  sweep  voltage  which  represents 
the  conventional  approach  to  Langmuir  probe  operation  wherein 
some  form  of  continuous  voltage  sweep  is  applied  between 
voltage  limits  V_  and  V+.  Fig.  1(B)  shows  the  pulse-modulated 

3 

sweep  which  has  been  utilized  with  P  .  The  voltage  pulses 
which  follow  the  sawtooth  envelope  generate  the  probe's 
current-voltage  characteristic.  During  the  interpulse 
period,  at  constant  voltage  V  ,  the  collected  probe  current 

D 

IB  provides  a  direct  measure  of  variations  in  the  probe- 
plasma  system.  The  pulse  duty-cycle  is  short  so  that  the 


2 


probe  rests  at  its  baseline  potential  Vg  for  a  period  much 
longer  than  the  pulse  width.  Fig.  2  shows  pulse  and  base¬ 
line  durations  identified  as  T  and  T„ ,  respectively.  So 

on  B 

that  sweep  voltage  transients  will  not  affect  the  value  of 
probe  current,  the  probe  current  is  sampled  at  the  termination 
of  subinterval  within  the  sweep  pulse  and  at  the  center 
of  the  baseline  interval  T_. 

D 

With  T  much  shorter  than  either  T_  or  the  time  constant 
on  B 

4 

of  the  probe  surface  contamination  layer  ,  the  pulse  procedure 
will  maintain  the  surface  condition  and  associated  voltage 
drop  at  a  more  nearly  constant  level  than  when  using  a 
continuous,  slowly-varying  sweep  voltage.  The  resulting 
current-voltage  characteristic  can  then  be  unfolded  from  the 
plasma  density  fluctuations  (5Ig)  so  that  the  electron 
temperature  and  density  are  determined  uniquely.  In  addition, 
the  I  -values  provide  the  raw  data  from  which  density 

O 

fluctuation  power  spectra  are  determined. 

3 

Payload  configuration  and  P  characteristics.  A  pair 
of  pulsed  probes  were  diametrically  extended  from  the 
forward  end  of  the  rocket  payload  (Fig.  3).  The  sensing 
elements  constructed  from  tungsten  wire,  were  isolated  from 
their  extension  booms  by  coaxial  guard  electrodes  driven  at 
the  same  potential  as  the  probes  themselves.  One  of  the 
probes,  defined  as  the  I-probe,  operated  with  Vg^-lv, 
yielding  net  ion  baseline  current  I„.  The  other  probe, 
defined  as  the  E-probe,  operated  with  V  n.  +  2v,  yielding 
net  electron  baseline  current  Ig.  Both  probes  generated 
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complete  current-voltage  characteristics  in  Tg  £  400  msec, 
yielding  absolute  values  of  and  Tg  at  an  approximate  2.5 
Hz  rate.  Maximum  Ig  sampling  occurred  at  2048  Hz,  resulting 
in  0.5  meter  spatial  resolution  for  relative  electron  density 
fluctuations  at  a  vehicle  velocity  of  1  km/sec. 

Probe  electrometers  were  set  to  operate  over  a  dynamic 

-10  -4 

range  extending  from  4(10  )  to  2.5  (10  )  amperes,  with 

automatic  switching  over  8  ranges  maintaining  9  bit  accuracy 
for  all  anticipated  ionospheric  conditions.  The  automatic 
ranging  is  best  illustrated  in  Figure  4,  an  actual  in-flight 
analog  record  of  telemetry  channel  outputs  for  the  probe 
currents  and  applied  voltages,  the  roll  magnetometer  for 
magnetic  aspect  determinations,  and  pitch,  yaw  and  roll 
monitors  on  the  ACS  jets.  The  data  sample  presented  in 
Figure  4  was  collected  when  the  payload  was  55  seconds  into 
flight  while  the  probe  electrometers  were  being  driven 
through  a  load  resistor  for  calibration.  The  probes' 
operation  alternated  between  a  fixed-bias  mode  and  a  pulsed- 
sweep  mode,  with  absolute  currents  determined  by  a  simple 
algorithm  which  coupled  the  switching  0-5  v  TM  signal  on  PCM 
channel  26-1  with  the  sweep  current  range  monitor  on  channel 
27-1.  The  record  format  in  Fig.  4  helped  provide  field- 
estimates  of  density  profiles  without  distortions  of  magnetic 
aspect  sensitivities  and  attitude  control  jets. 

III.  RESULTS 

Density  profiles  and  irregularity  structures.  By 
9  P.M.  (LT)  on  the  night  of  the  first  rocket  launch  (PLUMEX  I) 
the  bottomside  of  the  F-region  had  risen  to  an  approximate 
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altitude  of  400  km.  The  F-region  then  began  a  downward 
drift  with  a  simultaneous  onset  of  spread-F.  The  downward 
drifting  and  spread-F  conditions  continued,  and  when  the 
bottomside  F-layer  had  descended  to  an  altitude  below  300 
km,  the  rocket  was  launched  (12:31:30  UT  on  day  198;  00:31:30, 
17  July  1979,  LT) . 

Figure  3  displays  the  upleg  measurements  of  relative 
electron  density  as  presented  by  correlated  ion-  (1„)  and 
electron-  saturation  I^(V+)  currents.  The  ordinate  has  a 
linear  scale  for  t ime-a f ter- launc h  with  altitude  super¬ 
imposed  at  50  second  increments.  Because  ion  and  electron 

saturation  currents  have  significantly  different  sensitivities 

8 

to  velocity,  sheath  and  magnetic  field  effects  ,  variations 
i  i.  4* 

in  1^  and  lg(V  )  not  mutually  corroborated  were  attributed 
to  the  various  aspect  sensitivities  and  excluded  from  Figure 
5.  This  approach  facilitated  analysis,  reduced  computer 
time,  and  established  credibility  in  the  interpretation  of 
the  curves  as  relative  electron  density  profiles. 

The  results  in  Figure  5  show  that  a  number  of  major 
depletions  (ANg/N°  <  0.9)  were  distributed  throughout  the  F- 
region.  Each  of  the  large  scale  depletions  (identified 
alphabetically)  has  its  own  distribution  of  irregularities, 
illustrated  in  Fig.  6  by  the  expanded  view  of  regions  C,  D, 

H  and  I.  It  is  clear  that  "C"  is  not  a  single  narrow  bite- 
out  but  a  collection  of  rather  large  irregular  structures 
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extending  over  a  total  altitude  domain  of  about  12  k-  . 
(Vehicle  velocity  in  region  C  was  2.4  km/sec.)  To  develop  a 
quantitative  view  of  irregularity  fluctuations  observed  in 
the  F -region,  continuous  linear  detrends  were  executed 
throughout  the  entire  upleg  trajectory.  The  variations 
about  those  linear  detrends  were  then  plotted  in  Figure  5  as 
"Irregularity  Intensity",  with  a  maximum  relative  scale  of 
+4.  A  fluctuation  as  great  as  +  4  approximately  represents 
a  +  80%  fluctuation  about  the  linear  detrend.  (Correlation 
of  these  results  with  macroscale  gradients  and  Altair  back- 
scatter  contours  are  discussed  in  a  companion  paper^ . 

Absolute  density  and  temperature ■  Absolute  values  of 
electron  density  and  temperature  were  determined  by  conven- 

9 

tional  analyses  of  Langmuir  probe  characteristics  with 
appropriate  care  to  eliminate  perturbing  effects  of  surface 
contamination^,  density  fluctuations^’10  and  magnetic  field 

g 

effects  .  Ana 1 y s is  of  approximately  25  characteristics  were 

executed  over  the  F-layer  from  340-560  km.  In  each  case  a 

-  3  + 

conversion  coefficient  a  =  Ng  [cm  ]/l  (V  )  was  determined 

so  that  the  Ig(V  )  profile  in  Figure  5  could  be  directly 

scaled  to  absolute  electron  densities.  This  procedure 

yielded  a  =  (5.5  +  0.5)  101,11  electrons  cm  ^A  ^ . 

The  upleg  profile  has  been  reconstructed  in  Figure  7 

with  relative  and  absolute  electron  density  plotted  as  a 

function  of  altitude.  The  result  shows  the  F-peak  at  375 

6  —3 

km,  with  a  maximum  density  of  1.3  (10  )cm  (+  10%). 
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Analysis  of  the  re  tarding-f  ield  region  of  the  same  set 


of  current-voltage  characteristics  yielded  I  =  (1350  + 

2  5  0 )  °  K  ,  with  no  obvious  signatures  of  electron  energy 
redistribution  in  and  around  the  depletions. 

Intermediate  wavelength  power  spectra .  The  pulsed 
probe  data  provided  an  excellent  opportunity  for  comparison 
with  the  numerical  simulations^  of  the  collisional  Rayleigh- 
Taylor  (R-T)  instability  at  intermediate  wavelengths. 
Attention  is  focused  on  the  bottomside  F-layer  gradient  and 
region  C,  which  is  believed  representative  of  the  mid-phase 
development  of  the  R-T  process'^.  Typically,  computer 
simulations  employ  several  values  for  the  zero-order  gradient 
scale  length 


and  initialize  the  code  with  some  two-dimensional  perturbation 
superimposed.  In  the  work  of  Keskinen,  et  al . ^  L  was 

selected  at  5,  10  and  15  km  and  the  perturbation  took  the 

,  12 
form 


fiN  (x,y,t-0)  =  _A 

-  (1C  )  Sin  (k  y)  Cos  (k  x)  +  2(10  )  Sin  (2k  y) 

N°  y  x  y 

e 

with  k  and  k  being  the  horizontal  and  vertical  wavenumbers, 
x  y 

respectively.  Both  k^  and  k^  were  set  equal  to  2tt/960  m  in 
the  simulation.  In  addition,  the  computation  assumed  that  L 
was  centered  at  300  km. 
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Under  actual  conditions  encountered  in  PLUM K X  1  (Pig  7), 
the  bottoraside  F  -layer  gradient  extended  from  240  to  290  km. 

The  question  of  gradient  scale  length  can  be  studied  in 
Figure  8  where  it  is  shown  that  the  bottomside  gradient 
(encompassed  in  the  105-125  sec  time  frame)  is  not  characterized 
by  a  single  value  of  L.  In  region  "C"  (114s  <  t  <  122s)  L 

'Xj  f\j 

is  seen  to  vary  between  2  and  10  km,  whereas  adjacent  domains 
(110s  <  t  <  113s  and  122s  <  t  <  126s)  can  be  characterized 
by  L  =  25  km.  We  would  suggest  that  the  adjacent  domains 
are  representative  of  the  zero-order  gradient  scale  length 
and  that  L  =  25  km  would  be  a  more  appropriate  value  in  the 
numerical  simulation. 

In  any  event,  computer  simulations11  with  L  =  5,  10  and 

15  km  showed  that  linearly  unstable  modes  saturate  by  non¬ 
linear  generation  of  vertical  modes.  The  results  yield  one¬ 
dimensional  power  laws  (horizontal  and  vertical)  that  vary 

with  a  spectral  index  ( =  n  in  a  k  n)  between  2.0  and  2.5. 

e 

To  explore  this  result  within  the  context  of  region  "C", 

power  spectral  analyses  were  conducted  over  sliding  intervals 

of  2.4  km.  The  results,  presented  in  Figure  9,  show  that 

the  dominant  behavior  is  k  over  the  range  k  =  2iT/lkm  to 

—  1  8  5 

k  =  2TT/2.5m.  The  k  '  behavior  at  t  =  116.001  sec  is  a 
result  of  the  very  sharp  density  gradient  (see  region  "C" 

Fig.  6)  encompassed  by  the  domain  of  the  spectral  analysis. 

In  general  we  would  conclude  that  our  results  support 
the  numerical  simulations  of  Keskinen,  et  al.11.  We  do 


point  out  however  that  a  spectral  index  variation  from  2  to 
2.5  is  a  rather  broad  domain.  Further  testing  of  this 
support  can  be  achieved  with  an  L  =  25  km  simulation  and  a 
downward  drifting  F-layer  model  that  is  more  in  keeping  with 
the  actual  experimental  conditions.  The  F-layer  time- 
history  can  be  important  since  unstable  modes  appear  to 
require  times  in  excess  of  4,000  seconds  to  saturate.  ..a 
time  during  which  the  F-layer  encountered  in  PLUMEX  1 
drifted  downward  in  excess  of  40  km. 

IV.  SUMMARY 

In  each  of  the  two  rocket  operations  (PLUMEX  I  &  II) 
conducted  at  the  Kwajalein  Atoll  during  July  1979,  the  NRL 
pulsed  plasma  probe  performed  flawlessly.  General  results 
concerning  coordinated  rocket  and  radar  measurements  of 

small  and  large  scale  irregularities  have  been  discussed  in 

1  2 

companion  papers  *  .  Complementary  results  developed  here 

include : 

(a)  In  PLUMEX  I  a  number  of  major  depletions  (ANe/N°<  90%) 

were  distributed  throughout  the  F-region,  from  the  bottom- 

side  gradient  centered  near  260  km,  through  the  F-peak,  to  a 

topside  altitude  of  500  km.  The  most  intense  "in  situ" 

irregularities  occurred  on  the  bottomside  ledge  where  gradient 

scale  lengths  were  found  to  vary  between  2  and  25  km.  The 

power  spectral  density  in  this  region  of  intense  irregularities 

-2  5 

on  the  bottomside  was  dominated  by  a  k  power  law  over 

the  intermediate  wavelength  domain  k  =  2Tr/lkm  to  k  =  2Tr/2.5m. 

The  experimental  conditions  were  reasonably  matched  to  the 


numerical  simulations  of  Keskinen  et  al. 


,  and  a  comparison 


11 

of  the  two  resulted  in  general  agreement. 

max 

(b)  In  PLUMEX  I,  the  F-peak  was  at  375  km,  with  Ng 

1.3  ( 1 0 6 )  cm-3  (+  10%).  The  electron  energy  distribution 

was  characterized  by  T  =  (1350  +  250)°K  with  no  obvious 

e  - 

signatures  of  energy  redistribution  in  and  around  the 
deplet ions . 
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Fig.  3  —  PLUMEX  payload  configuration  and  nominal  trajectory. 
(This  figure  has  been  adapted  from  a  Sandia  report.) 
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Fig.  5  —  Relative  electron  density  profile  of  macroscale  features  as  measured  simultaneously 
by  ion  and  electron  saturation  probe  currents  collected  on  the  upleg  trajectory  of  PLUMEX  I. 
The  “irregularity  intensity”  provides  a  measure  of  smaller  scale  structure  with  a  ±  4  intensity 
approximately  equal  to  ±  80%  fluctuations  about  a  linear  detrend. 
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CAMBRIDGE,  MA  02139 

OICY  ATTN  D.B.  COX 
OICY  ATTN  J.P.  GILMORE 

COMPUTER  SCIENCES  CORPORATION 
6565  ARLINGTON  BLVD 
FALLS  CHURCH,  VA  22046 
OICY  ATTN  H.  BLANK 
OICY  ATTN  JOHN  SPOOR 
OICY  ATTN  C.  NAIL 

COMSAT  LABORATORIES 
LINTHICUM  ROAD 
CLARKSBURG,  MD  20734 

OICY  ATTN  G.  HYDE 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 

OICY  ATTN  H.  LOGSTON 

OICY  ATTN  SECURITY  (PAUL  PHILLIPS) 

ESL  INCL. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 

OICY  ATTN  J.  ROBERTS 
OICY  ATTN  JAMES  MARSHALL 
OICY  ATTN  C.W.  PRETTIE 
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FORD  AEROSPACE  AND  COMMUNICATIONS  CORP 
3939  FABIAN  WAY 
PALO  ALOT,  CA  94303 

0 ICY  ATTN  J.T.  MATTINGLEY 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 

0 ICY  ATTN  M.H.  BORTNER  SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.O.  BOX  1122 
SYRACUSE,  NY  13201 

0  ICY  ATTN  F.  RE  I  BE  RT 

GENERAL  ELECTRIC  COMPANY 
TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O.  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
OICY  ATTN  DAS I AC 
01CY  ATTN  DON  CHANDLER 
OICY  ATTN  TOM  BARRETT 
OICY  ATTN  TIM  STEPHANS 
OICY  ATTN  WARREN  S.  KNAPP 
OICY  ATTN  WILLIAM  MCNAMARA 
OICY  ATTN  B.  GAMBIIL 
OICY  ATTN  MACK  STANTON 

GENERAL  ELECTRIC  TECH  SERVICES  CO.,  INC. 

HMES 

COURT 

SYRACUSE,  NY  13201 

OICY  ATTN  G.  MILLMAN 

GENERAL  RESEARCH  CORPORATION 
SANTA  BARBARA  DIVISION 
P.O.  BOX  6770 
SANTA  BARBARA,  CA  93111 

OICY  ATTN  JOE  ISE  JR 
OICY  ATTN  JOEL  GARBARINO 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  A K  99  701 

(ALL  CLASS  ATTN:  SECURITY  OFFICERS) 
OICY  ATTN  T.N.  DAVIS  (UNCL  ONLY) 

OICY  ATTN  NEAL  BROWN  (UNCL  ONLY) 

OICY  ATTN  TECHNICAL  LIBRARY 


GTE  SYLVANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

OICY  ATTN  MARSHAL  CROSS 

ILLINOIS,  UNIVERSITY  OF 
DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
URBANA,  IL  61803 

OICY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
400  ARMY -NAVY  DRIVE 
ARLINGTON,  V A  22202 

OICY  ATTN  J.M.  A£ IN 
OICY  ATTN  HANS  WOLFHARD 
OICY  ATTN  JOEL  BENGSTON 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

OICY  ATTN  DONALD  HANSEN 

INTL  TEL  &  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

OICY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

1401  CAMINO  DEL  MAR 
DEL  MAR,  CA  92014 

OICY  ATTN  S.R.  GOLDMAN 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

OICY  ATTN  DOCUMENT  LIBRARIAN 
OICY  ATTN  THOMAS  POTEMRA 
OICY  ATTN  JOHN  DASSOULAS 

LOCKHEED  MISSILES  &  SPACE  CO  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 

OICY  ATTN  DEPT  60-12 
OICY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  AND  SPACE  CO  INC 
3251  HANOVER  STREET 
PALO  ALOT,  CA  94304 

OICY  ATTN  MARTIN  WALT  DEPT  52-10 
OICY  ATTN  RICHARD  G.  JOHNSON  DEPT 
52-12 

OICY  ATTN  W.L.  IMHOF  DEPT  52-12 


KAMAN  SCIENCES  COR? 

P.O.  BOX  7463 

COLORADO  SPRINGS ,  CO  30933 
OICY  ATTN  T.  MEAGHER 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 

0 ICY  ATTN  IRWIN  JACOBS 
0  ICY  AITN  I.  ROTHMUELLER 

LOWELL  RSCH  FOUNDATION,  UNIVERSITY  OF 
450  AIKEN  STREET 
LOWELL,  MA  01854 

0 ICY  ATTN  K.  BIBL 


M.I.T.  LINCOLN 
P.O.  BOX  73 

LABORATORY 

LEXINGTON, 

MA 

02173 

OICY 

ATTN 

DAVID  M.  TOWLE 

OICY 

ATTN 

p.  waldron 

OICY 

ATTN 

L.  LOUGHLIN 

OICY 

ATTN 

D.  CLARK 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  58  37 
ORLANDO,  FL  32805 

0 ICY  ATTN  R.  HEFFNER 

MCDONNELL  DOUGLAS  CORPORATION 
5  301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
0  ICY  ATTN  N.  HARRIS 
01CY  ATTN  J.  MOULE 
0 ICY  ATTN  GEORGE  MROZ 
01CY  ATTN  W.  OLSON 
01CY  ATTN  R.W.  HALPRIN 
01CY  ATTN  TECHNICAL  LIBRARY  SERVICES 

MISSION  RESEARCH  CORPORATION 

735  STATE  STREET 

SANTA  BARBARA,  CA  93101 

OICY  ATTN  P.  FISCHER 

01CY  ATTN  W.F.  CREVIER 

OICY  ATTN  STEVEN  L.  GUTSCHE 

OICY  ATTN  D.  SAPPENFIELD 

OICY  ATTN  R.  BOGUSCH 

OICY  ATTN  RALPH  KILB 

OICY  ATTN  R.  HENDRICK 

OICY  ATTN  DAVE  SOWLE 

OICY  ATTN  F.  FAJEN 

OICY  ATTN  M.  SCHEIBE 

OICY  ATTN  CONRAD  L.  LONGMIRE 

OICY  ATTN  WARREN  A.  SCHLUETER 


MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 

OICY  ATTN  JOHN  MORGANSTERN 
OICY  ATTN  G.  HARDING 
OICY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 

OICY  ATTN  W.  HALL 
OICY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
1456  CLOVE RFIELD  BLVD. 

SANTA  MONICA,  CA  90404 

OICY  ATTN  E.C.  FIELD  JR 

PENNSYLVANIA  STATE  UNIVERSITY 
IONISPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 

(NO  CLASSIFIED  TO  THIS  ADDRESS) 

OICY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  INC. 

442  MARRETT  ROAD 
LEXINGTON,  MA  02173 

OICY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 

OICY  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS  INC. 

P.O.  BOX  1069 
BERKELEY,  CA  94701 

OICY  ATTN  A.  THOMPSON 

R  &  D  ASSOCIATES 

P.O.  BOX  9695 

MARINA  DEL  REY,  CA  90291 

OICY  ATTN  FORREST  GILMORE 
OICY  ATTN  BRYAN  GABBARD 
OICY  ATTN  WILLIAM  B.  WRIGHT  JR. 
OICY  ATTN  WILLIAM  J.  KARZAS 
OICY  ATTN  ROBERT  F.  LELEVIER 
OICY  ATTN  H.  ORY 
OICY  ATTN  C.  MACDONALD 
OICY  ATTN  R.  TURCO 


RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 

OICY  ATTN  CULLEN  CRAIN 
0 ICY  ATTN  ED  BEDROZIAN 

RIVERSIDE  RESEARCH  INSTITUTE 
80  WEST  END  AVENUE 
NEW  YORK,  NY  10023 

01CY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATION,  INC. 

P.O.  BOX  2351 
LAJOLLA,  CA  92038 

01CY  ATTN  LEWIS  M.  LINSON 
0 ICY  ATTN  DANIEL  A.  HAMLIN 
01CY  ATTN  D.  SACHS 
01CY  ATTN  E.A.  STRAKER 
01CY  ATTN  CURTUS  A.  SMITH 
OICY  ATTN  JACK  MCDOUGALL 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

OICY  ATTN  BARBARA  ADAMS 

SCIENCE  APPLICATIONS,  INC. 
HUNTSVILLE  DIVISION 
2109  W.  CLINTON  AVENUE 
SUITE  700 

HUNTSVILLE,  AL  35805 

OICY  ATTN  DALE  H.  DAVIS 

SCIENCE  APPLICATIONS,  INCORPORATED 
8400  WESTPARK  DRIVE 
MCLEAN,  VA  22101 

OICY  ATTN  J.  COCKAYNE 

SCIENCE  APPLICATIONS,  INC. 

80  MISSION  DRIVE 
PLEASANTON,  CA  94566 
OICY  ATTN  SZ 

SRI  INTERNATIONAL 
3  33  RAVE  NS  WOOD  AVENUE 
MENLO  PARK,  CA  94025 

OICY  ATTN  DONARD  NEILSON 
OICY  ATTN  ALAN  BURNS 
OICY  ATTN  G.  SMITH 
OICY  ATTN  L.L.  COBB 
OICY  DAVID  A.  JOHNSON 
OICY  ATTN  WALTER  G.  CHESNUT 
OICY  ATTN  CHARLES  L.  RINO 
OICY  ATTN  WALTER  JAYE 
OICY  ATTN  M.  BARON 


OICY  ATTN  RAY  L.  LEADABRAND 
OICY  ATTN  G.  CARPENTER 
OICY  ATTN  G.  PRICE 
OICY  ATTN  J.  PETERSON 
OICY  ATTN  R.  HAKE,  JR. 

OICY  ATTN  V.  GONZALES 
OICY  ATTN  D.  MCDANIEL 
OICY  ATTN  R.  TSUNODA 
TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  017  30 

OICY  ATTN  W.P.  BOQUIST 

TRW  DEFENSE  &  SPACE  SYS  GROUP 

ONE  SPACE  PARK 

REDONDO  BEACH,  CA  90278 

OICY  ATTN  R.K.  PLEBUCH 
OICY  ATTN  S.  ALTSCHULER 
OICY  ATTN  D.  DEE 

VISIDYNE,  INC. 

19  THIRD  AVENUE 

NORTH  WEST  INDUSTRIAL  PARK 

BURLINGTON,  MA  01802 

OICY  ATTN  CHARLES  HUMPHREY 
OICY  ATTN  J.W.  CARPENTER 
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.  14850 
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DR. 

W. 

E.  SWARTZ 
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R. 
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R. 

HARTEL./CODE  621 

Washington,  DC  20350 
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INSTITUTE  FOR  DEFENSE  ANALYSIS 
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